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G^)pffipn  = Fffipn  /^Hn  for  Hr^/2)^-H^Hm/2) 

or 

G(CJ)p_p„  = 0 for  iO  otherwise. 

. mn 

snd  r„,  are  evaluated  at  U)  = Cd  in  Eq.  (4.4).  The 
vni  Wlw  m o 

n = m terms  are  the  covariance  of  the  same  mode  observed  at 
the  two  detectors.  With  the  range  of  tine  delay  |7^  » { 

limited  to  about  the  travel  time  between  the  detectors 
|(x1  - %2 )/Gm I » the  most  important  terms  in  the  summation  are 
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ABSTRACT 


The  usual  way  to  determine  the  direction  of  a radiat 
ing  source  using  an  array  is  to  steer  the  array  for  maximum 
output.  The  steering  can  be  done  by  time  delays  or  by  mechani 
cally  turning  the  array.  The  output  of  the  array  can  be  ex- 
pressed as  the  sum  of  the  time  average  of  the  products  of  the 
pressure  observed  at  detectors  m and  n . The  value 


of  <^ptt<Tm)  can  be  maximized  by  a proper  choice  of  time 

delays  T^»  ° This  procedure  is  straightforward  in  an  infi- 

nite homogeneous  medium.  If  the  medium  Is  a layered  waveguide, 
there  are  many  more  possibilities  for  submaxima  of  the  P_B^ 


The  normal-mode  solution  of  the  radiation  field  of  a 


band  limited  noise  point  source  in  a layered  waveguide  was 
given  by  the  author  in  J.  Acoust.  Soc.  Am.  31»  1^73-^79  0 959) 
The  value  of  r„cr„v  as  a function  of  time  delay  is 


compared  with  the  value  of  ?nPB  obtained  with  mechanical 
steering.  For  these  calculations  the  noise  source  is  assumed 
to  have  a bandwidth  of  1/15  of  the  center  frequency,  and  the 


depth  of  the  water  is  assumed  to  be  about  2^.  over  a thick  layer 
of  unconsolidated  sediment.  The  number  of  maxima  of  the 


Pn^nv  is  re*Late^  to  the  number  of  modes  propagating  in 
the  waveguide  if  the  steering  is  done  with  time  delays.  Meehan 
leal  steering  yields  one  maximum  that  corresponds  to  the  \source 
direction.  \ 


a 


1 « Introduction 

The  use  of  an  array  to  determine  the  direction  of  a 

radiating  source  has  been  studied  with  many  array  designs# 

(1  -7) 

processing  techniques,  end  signal-to-noise  conditions. 

The  usual  way  to  determine  the  direction  of  the  source  is  to 
adjust  the  phases  of  the  signals  observed  at  each  detector  so 
that  the  time  average  of  the  square  of  the  sums  of  the  outputs 
is  a maximum.  The  phase  adjustments  can  be  made  by  mechanically 
turning  the  array  or  by  time  delays.  Assumptions  common  to  all 
the  studies  are  that  the  array  is  in  a homogeneous  medium,  and, 
in  the  absence  of  noise,  that  the  signal  observed  at  one  detec- 
tor is  essentially  the  same  as  the  signal  observed  at  any  other 
detector.  If  the  medium  is  not  homogeneous  but  is  a layered 
waveguide,  a number  cf  things  can  happen.  The  medium  is  disper- 
sive, i.e0,  different  frequencies  travel  with  different  group 
and  phase  velocities.  Depending  on  the  number  of  solutions  al- 
lowed by  the  boundary  conditions,  a particular  frequency  can 
have  several  different  group  and  phase  velocities.  In  addition 
to  these  difficulties,  the  phase  and  group  velocities  are  not 
equal.  Thus  the  signal  changes  as  it  travels  in  the  guide,  and 
characteristics  of  the  received  signal  are  largely  due  to  the 
waveguide. 

In  this  paper  we  re-examine  the  array  steering  problem 
with  both  the  array  and  source  in  a layered  waveguide.  The  so- 


lution of  the  pressure  field  due  to  a point  source  that  radiates 
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band  limited  noise  in  a layered  waveguide  has  been  given  in 
terms  of  normal  modes.  This  pressure  field  is  used  with 

the  usual  procedure  of  array  calculations  in  the  following  sec- 
tions. A numerical  example  with  a two-element  array  and  a 
shallow-water  waveguide  is  given.  I should  remark  here  that 
this  development  is  mainly  concerned  with  waveguides  that  are 
a few  acoustic  wavelengths  thick  and  with  very  large  source- 
receiver  distances.  The  source  is  driven  by  band  limited  noise „ 
and,  in  the  numerical  example,  the  bandwidth  is  one  fifteenth 
of  the  center  frequency. 

2.  Propagation  in  a Layered  Waveguide 

Tolstoy  has  given  the  theory  of  propagation  of  acous- 
tic waves  due  to  a simple  harmonic  point  source  in  a layered 
waveguide  in  detail. ls  sufficient  here  to  summarize 
his  results. 

All  energy  from  a source  in  a uniform  layer  that  is 
totally  reflected  at  the  boundaries  of  the  layer  is  trapped  in 
the  layer.  The  trapped  energy  spreads  cylindrically  and  has 
many  reflections  from  the  top  of,  within,  and  at  the  bottom  of 
the  waveguide.  The  most  convenient  way  to  describe  the  pressure 
field  at  a large  distance  is  to  use  the  normal-mode  formulation 
of  the  problem.  We  assume  a horizontally  stratified  medium  with 
a free  surface  on  the  top.  This  is  basically  the  shallow-water 
problem  with  water  over  layers  of  sediment  extending  to  basement. 

The  simple  harmonic  source  is  at  depth  d in  the 
first  layer.  The  detector  is  at  depth  z in  the  first  layer 
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and  at  range  R from  the  source.  The  long-range  solution  for 
the  acoustical  pressure  P;  Is  the  following: 

(2.1  ) 

M 

P;  = iPlAR“1/2  e^ut  Z pm  sin  r1ffl  z sin  r1m  d x exp  f-i  (KmF+ff> ] 

m— I 


in  which 


CO  - the  angular  frequency , 

Km  ~ the  horizontal  component  of  the  wave  number  for 
the  mth  mode, 

K = w/c  , 

Cm  = the  phase  velocity  of  the  mth  mode, 
rim  ~ vei>tical  component  of  the  wave  number  in  the 
1th  layer, 

rim  " ^2/“l2  • 

Pm  55  the  mode  excitation  for  mth  mode, 

= velocity  of  sound  in  the  ith  layer. 


pA  = density  in  the  1th  layer. 

It  Is  difficult  to  express  /(  . r._,  and  p_  as  func- 

ra  xin  in 

tlons  of  CJ  in  terms  of  the  parameters  of  the  waveguide.  For 
a particular  waveguide,  numerical  solutions  are  computed.  There 
are  many  papers  in  the  literature  that  pive  the  theory  and  tech- 
niques. 1-1  6 5 For  our  purpose,  the  most  important  thing  is 
that  simple  functions  of  UJ  can  be  used  to  approximate  the  nu- 
merical values  of  rlm(ui),  Km(w),  and  pffl<io)  . 
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The  extension  of  the  simple  harmonic  source  to  the 
band  Diraited  noise  source  Is  obtained  from  Fq.  (2.1)  by  using 
the  Fourier  Integral  expression  for  the  noise  source: 

(2.2) 

P(t)-[  /g(6o)lAB-1^2e1^t  Z pmsin  r1mz  sin  r1rad  * expf * dtd  • 
-co 1 111=1  J 

The  source  and  receiver  filter  characteristics  are  Included  In 

the  function  g(cO)  . 

3.  Array  of  Detectors 

The  usual  vay  to  process  the  signals  from  an  array 
of  detectors  is  to  add  the  signals  and  take  the  time  average 
of  the  square  of  the  sum.  The  direction  to  source  can  be  de- 
termined by  measuring  the  output  of  the  array  as  a function  of 
steering  direction.  The  objective  of  steering  is  to  match  the 
phases  of  the  signal  from  each  detector  so  that  all  the  signals 
are  the  same.  Array  steering  can  be  done  by  turning  the  array 
or  by  placing  the  proper  time  delay  in  each  signal  channel. 

Let  us  assume  that  we  have  an  array  of  detectors  1, 

2,  ...  N . The  detectors  measure  acoustical  pressures  (t), 
P2(t),  ...  Pjj(t)  . Time  delays  are  placed  in  each  signal  chan- 
nel so  the  signals  observed  become  (without  change  of  notation) 


pi  (t-T^  ) $ P2(t-T2)*  •••  • The  outPUt  of  an  additive 
array  is  the  sum  of  the  individual  signals.  The  time  average 
of  the  square  of  the  array  output  is  the  following: 


\ / 


(3.1 ) 


N K . * v 

F = f f < • Pr>  (t-rn)  > 


in  which  P is  the  complex  conjugate  of  P and 


(3>2) 


( PmPn  ) - limT->  o©2T  pm(t“^m5Pn  (t"Vdt 


Fquation  (3«1 ) can  be  separated  into  two  sets  of  terms,  the 

squared  terms  ^ ^tnPra  ^ and  cross  ^<3rtT’s  ^ PmPn  ^ • 

The  squared  terms  are  independent  of  steering  and  time  delay 
for  a source  of  stationary  noise.  The  directional  information 
is  obtained  from  cross  terms.  Upon  separation  of  the  term*  in 
the  summation,  Fq<  (3.1 ) becomes 


(3.3) 


in  which 
(3.M 


N N 


F "f  < p„2>  + f f VVfn’  • 

n/m 

<W : <V*-V  • F„*<t-rn)) 


The  cross  terras  are  also  the  covariances  of  pairs  of 
signals  with  time  delays  as  defined  by  Eq.  (3«4).  The  bal- 

ance of  this  paper  will  discuss  a single  covariance  F.12  . This 
simplification  is  mainly  one  of  notation.  It  does  not  limit  the 
results  to  a two-element  array  because  all  of  the  cross  terms 
in  Eq.  (3.3)  are  similar. 

4.  Covariance  of  Two  Signals 
The  coordinate  axis  is  shown  in  Fig.  1.  The  detectors 
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are  in  the  xz  plane  at  (x^  ,z^  ) and  ^2*^2)  0 source 

is  at  depth  d * distance  3 » and  angle  0 as  measured  from 
the  y axis.  The  distance  of  the  source  is  presumed  to  be 


large  relative  to  the  thickness  of  the  waveguide  and  the  sepa- 
ration between  the  detectors  * 

The  covariance  is  obtained  from  Fq,  (2„2)  and  Eq.  (3.2 )„ 
The  product  of  the  integrals  can  be  transformed  to  £ multiple 
integral  by  changing  the  variable  of  integration  fron  CJ  to 
L)1  . The  integration  over  t and  application  of  the  limit  as 


T tends  to  infinity  indicates  that  the  expression  vanishes  for 
U/u'  Thus  one  obtains  the  followings 


(4.1  ) 


CO 

/ cU))Pn,Pnsin  r,nzl5ln  r,nz2 


in  which 


x sin  r1md  sin  r1nd  exp  ( ~'r2 ) ]}  dW  » 

o(cj)  = g g 

R1  5 R2~2Rx.j  sin  e + x- 2 
Rg  = R2-2Rx2  sin  0 + x?2 
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* 


L 


(4.2) 


Rj  ~ R -x1  sin  e 


E » x. 


R » x. 


Rg  ~ H -Xg  sin  e 


Equation  (4.1)  can  be  evaluated  if  the  dependences 


of  G,  P , r , and  K.  upon  6 J are  known.  Numerical  values 
Dl  10  m 


of  K and  r as  functions  of  CJ  are  shown  in  Fig.  2A„ 
in  id 


and  numerical  values  of  pm(cj)  are  shown  in  Fig.  2B.  The 


calculations  were  made  by  Tolstoy  for  the  shallow-water  section 

174  (18) 


near  Fire  Island.  J For  first  approximation*  /Cm  and  rffi 
may  be  expanded  as  linear  functions  of  0J  . Using  the  follow- 


ing definition  of  the  group  velocity 


= <&L 


« d/cn 


we  have 


(4.3) 


K.  ^ 


(W~a  )/U 
^ nr'  m 


r»  * v-b. 


for  (j)  > 0 


The  approximate  expressions  for  and  rffl  can  be 

substituted  into  Eq.  (4.1).  If  the  bandwidth  of  the  filter 


function  G(u>)  is  narrow,  a simple  function  can  be  used  for 


p and  the  integrals  evaluated.  Wide  bands  can  be  treated  by 
m 


piecewise  integration  of  several  narrow  bands.  One  must  be 


careful  in  this  integration  since  and  K are  positive 


for  outgoing  waves,  (t)  is  restricted  to  positive  values  in 


,lk  ^ — 17 


IN  CPS 


Fq.  (4.3). 

The  bandwidth  of  interest  is  the  product  of  the  filter 
function  G (w)  and  the  mode  excitations  PraPn  • The  product 
may  be  approximated  by  an  ideal  filter  with  a bandwidth  LCJ 
and  center  frequency  C0Q  . The  terms  in  F.q.  (4.1  ) can  be  in- 
tegrated with  Fqs.  (4.2)  and  (4.3).  The  change  of  rK  over 
the  range  of  is  small  enough  that  rffl  can  be  considered 

constant  for  this  example.  For  AW  about  ^/^5  (jjQ  • the  er- 
rors are  less  than  10  percent  of  the  mode  interference  terms. 


The  result  of  Integration  over  the  bandwidth  AU^n  is  the  fol= 


lowing : 
(4.4) 


USt T72  m?n  PmPn  sln  r1mz1  sin  r1nz2  sin  r1md  sin  r1nd 


R1  ' *2 


in  which 


slne'  ' 
0 = [«„  - /<„>  R - -Kn*2>  sin  e -U)a  j 


TT  - TT 
m n 


x eos  a 5ln 

008  P *0^0/2 


*1  *2 
JT  - TT 

m n 


w*w0  , 


nol-  = equivalent  ideal  filter  bendvldth  for  the  product 


0<0»pnp„  , 


Pn  = the  average  value  of  pm  in  the  frequency  band 
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and 


G<P)pmP„  = PmP„  //Hn  for 


or 

G((J)pmp^  = 0 for  otherwise, 
mn 

£L»  and  r are  evaluated  at  CO  - (a)  in  Eq.  (4.4).  The 
vui  m m o ^ 

n = m terms  are  the  covariance  of  the  same  mode  observed  at 
the  two  detectors.  With  the  range  of  time  delay  |Tj  - T,, 1 
limited  to  about  the  travel  time  between  the  detectors 
|(x1  - x2VVm\  9 the  most  imP°rtant  terms  in  the  summation  are 
the  n = m terms  for  very  large  source  distance.  The  cross  terms 
can  have  a contribution  for  time  delays  that  are  about  the  dif- 
ference of  the  travel  times  of  energy  in  the  different  modes. 

The  time  delays  that  maximize  the  covariance  of  the  cross  terms 
or  mode  interference  terms  are  related  to  the  distance  from  the 
source  to  the  receiver.  For  a single  detector,  the  auto- 
covariance  function  can  be  used  to  determine  the  source  distance. 

We  assume  in  the  array  steering  problem  that  the 
source  is  at  great  distance  and  that  the  time  delay  range  is 
limited  to  that  required  to  steer  the  array.  The  cross  terms, 
n ^ m , will  be  dropped.  With  these  assumptions  and  simplifi- 
cations Eq.  (4.4)  becomes 


(4.5) 


e12~  T6"  Z V sln  r1m*l  sln  r1mz2  rlB  d 


2 2 


X cos (wo'T-/Cmx  sin  6)  — 


3 in  ^r-  ) 


~^IT 


x sin  e 
Um 


— 1 2- 


in  which 

x1  = x/2  , 

Xg  =-x/2  , 

r = r2.r,  . 

R-p1/2R  1/2 

n — Kg  * 

and 

= equivalent  ideal  filter  bandwidth 

p 

for  the  product  G(6J)p_  . 

Ul 

The  mth  term  in  the  summation  is  the  covariance  of 

the  acoustical  pressures  at  the  detectors  for  the  mth  mode. 

The  size  of  each  term  is  dependent  upon  the  detector  positions 

and  source  depth,  in  addition  to  the  array  steering  parameters. 

The  maximum  value  of  the  contribution  of  the  mth  mode 

to  the  total  covariance  F,^  is  determined  by  the  mode  excita- 

— 2 2 

tion  pm  and  source  depth  d , i.e.,  the  value  of  pffl  sin  r1nja 
The  detector  depths  enter  as  the  product 
sin  r1p)  z1  sin  r1ffl  z2  . The  value  of  this  product  can  be  plus, 
minus,  or  zero.  Since  we  are  interested  in  the  covariance  of 
the  same  mode  at  two  detectors  and  the  total  covariance  as  a 
function  of  steering,  we  assume  that  the  detectors  are  at  the 
same  depth. 

The  dependence  of  the  envelope  of  the  covariance  func- 
tion upon  steering  is  given  by  ttK  term  of  form  sin£/£  . The 
maximum  of  sin£/£  is  at  £ s 0 . The  value  of  time  delay  X 
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for  which  the  function  is  maximum  is  dependent  upon  the  proup 
velocity  Uffl  (for  x sin  6 not  zero).  The  sin£/£  term 
multiplies  the  term  cosCcj sin  e)  that  is  elso  dependent 
upon  the  steering  and  mode.  This  term  gives  the  oscillation  of 
the  covariance  of  the  mth  mode  as  a function  of  steering.  It  is 
dependent  upon  the  angular  frequency  and  the  horizontal 

component  of  the  wave  number  K „ (or  in  terms  of  the  phase 

ID 

velocity,  Cm  = /<m/k>0  )•  The  phase  and  group  velocities  of 
energy  propagation  in  a waveguide  are  generally  different,  and 
even  for  a single  mode  the  steering  time  delay  T Tor  maximum 
of  the  covariance  may  not  correspond  to  the  direction  to  the 
source.  Since  the  group  velocity  U and  phase  velocity  C 

u.  ui 

are  dependent  upon  the  mode,  the  covariance  of  the  acoustical 
pressures  for  each  mode  is  a maximum  at  different  time  delays 
(for  x sin  6 not  zero).  If  T is  zero  and  sin  e is  varied, 
i.e.,  if  the  source  is  moved  across  the  acoustic  axis  of  the 
array  or  the  array  is  turned  for  a stationary  source,  then  all 
of  the  terms  are  maximum  at  x sin  6=0.  For  either  type  of 
steering,  the  sum  of  the  terms  is  complicated,  and  it  is  more 
instructive  to  consider  specific  numerical  examples. 

5.  Numerical  Example 

The  shsllow-water  area  off  Fire  Island  has  been  used 
for  studies  of  propagation  of  acoustic  waves  in  a layered  wave- 
guide. The  agreement  of  the  experimental  data  with  theoretical 
calculations  has  been  good.^®*  The  waveguide  layering  and 


the  curves  of  uJ  vs  the  horizontal  and  vertical  components  of 
wave  number  for  the  Fire  Island  area  are  shown  in  Figs.  1 and  2. 

The  source  is  assumed  to  be  about  7 ® deep.  The  de- 
tector3  are  assumed  to  be  500  m apart  and  19  ® deep.  The  center 
angular  frequency  of  the  source  Is  assumed  to  be  300ir  rad/sec 
with  bandwidth  of  20*  rad/sec  . With  these  numbers,  the  con- 
tribution to  the  covariance  of  the  1st,  2nd,  and  4th  modes  are 
nearly  equal,  and  the  3rd  mode  is  nearly  zero.  Numerical  values 
of  K _ ard  U are  obtained  from  Fig.  2.  The  amplitude  factors 
of  the  modes  are  nearly  equal.  The  normalized  covariance 
is  the  following: 

^ sin(zvy2)(r~x  sin  e/U ) 

F12<r>  = 1/3  ^,r.2  ^os<wor-V  sln  e)  W sin  e/8B) ' 


The  covariance  function  for  an  array  In  an  infinite  homogeneous 
medium  is  obtained  by  letting  the  phase  and  group  velocities  be 
equal  lor  all  modes.  For  the  last  case,  the  function  is  of  the 


following  form: 


where 


cos  cj 


sln  QQS/2 

*5? 


£ = U T - x sln  S/a1  i 


the  covariance  is  the  same  function  for  either  time  delay  or 
sin  6 steering. 

The  results  of  the  calculations  are  shown  in  Fig.  3- 
Fipure  3A  shows  the  covariance  of  two  elements  in  an  infinite 
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hcaofeneous  ocean.  This  result  is  similar  in  appearance  to  the 
'©variance  of  a single  mode  in  the  shallow-water  waveguide  for 
T = 0 and  sin  P steering.  Figure  3B  shows  the  total  covari- 
ance function  for  sin  6 steering  and  three  modes  propagating. 
The  sharp  maximum  at  x sin  6=0  and  the  lack  of  the  sin£/£ 
envelope  of  the  total  covariance  function  are  due  to  the  dif- 
ference between  the  K* s of  the  different  modes. 

m 

Figure  3C  shows  the  total  covariance  for  the  second 
numerical  example.  The  source  is  assumed  to  be  in  the  direction 
e = 90  deg  or  in  line  with  the  detectors.  The  covariance  func- 
tion for  each  of  the  modes  is  similar  in  appearance  to  that  in 
Fig.  3A.  However,  the  frequency  of  oscillation  is  different 
for  eech  of  the  modes,  and  the  phase  shifts  of  the  cosine  waves 
are  also  different.  The  result  of  this  is  that  the  covariance 
functions  for  the  modes  constructively  and  destructively  inter= 
fere  with  each  other.  The  phase  and  group  velocities  are  not 
equal,  and  the  phase  of  the  cosine  wave  is  not  zero  at  a steer- 
ing delay  such  that  the  sin^/f  term  is  maximum.  Since  the 
group  velocities  of  different  modes  are  different,  the  enve- 
lopes of  the  covariance  functions  sin  are  maximum  at 

different  values  of  time  delay.  The  positions  of  the  maxima 
of  the  envelopes  for  the  three  modes  are  indicated.  The  exact 
shape  of  the  interference  pattern  in  Fig.  3C  Is  very  much  de- 
pendent upon  the  source  direction  and  the  detpctor  separation 
x sin  0 . 
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If  we  were  to  compute  the  total  covariance  for  a small 


ranpe  of  e near  eQ  , we  would  find  that  the  interference  max- 
ima and  uiLnima  would  describe  an  envelope  that  is  approximately 
of  the  form  sin £/£  . This  suggests  that  the  total  covariance 
function  could  be  squared  and  averaged  over  o small  range  of 
x sin  e . (This  type  of  data  processing  does  not  seem  to  be 
advantageous  because  it  requires  many  determinations  of  the 
covariance  squared  as  x sin  e is  varied.)  It  may  be  possible 
to  combine  the  outputs  of  several  detectors  at  different  depths 
so  as  to  detect  each  mode  separately.  The  separate  covariance 
function  for  each  mode  could  be  examined  separately  or  recom- 
bined with  suitable  time  delays  to  give  a total  covariance 
similar  to  that  in  Fig.  3B. 

6.  Conclusions 

Steering  an  array  in  a waveguide  is  more  difficult 
than  steering  an  array  in  an  infinite  homogeneous  medium.  The 
reason  is  that  the  wave  propagation  in  a waveguide  is  more  com- 
plex. The  energy  propagates  in  many  modes,  and  each  mode  can 
have  different  phase  and  group  velocities.  Although  we  have 
not  studied  the  signal-to-incoherent-noise  problem,  it  is  rea- 
sonable to  compare  signal-to-noise  studies  with  the  array 
steering  problem.  In  both  cases  we  assume  that  wo  are  measuring 
the  covariance  of  the  acoustical  pressure  observed  at  two  de- 
tectors. A convenient  definition  of  the  output  signal-to-noiso 


h 


ratio  of  a covariance  detector  is  the  ratio  of  the  covariance 

(2 ) 

squared  to  the  wean  square  fluctuations  of  the  covariance  <> 

In  this  respect  then,  fluctuations  of  the  covariance  function 
for  very  small  changes  of  source  direction  can  be  considered 
as  an  apparent  noise. 

Mechanical  or  sin  0 steering  gives  a covariance  of 
the  acoustical  pressures  with  a well-defined  maximum.  The  co- 
variance  maximum  is  in  the  proper  direction.  The  width  of  the 
maximum  of  the  envelope  of  the  covariance  is  narrower  than 
would  be  expected  on  the  source  bandwidth.  This  might  be  ex- 
pected since  several  independent  pieces  of  information,  i.e., 
several  propagation  modes,  are  combined  in  this  calculation,, 

The  maximum  of  the  covariance  function  should  not  fluctuate 
for  small  changes  of  the  source  direction.  In  the  absence  of 
incoherent  background  noise,  the  apparent  noise  level  for  this 
type  of  array  steering  should  be  very  low. 

Time  delay  steering  of  a two-element  array  that  de- 
tects all  modes  yields  a covariance  that  has  several  maxima 
and  minima.  The  complication  is  caused  by  the  constructive 
and  destructive  interference  of  the  covariance  functions  of 
the  several  modes.  The  magnitude  of  the  interference  maxima 
and  minima  is  less  than  or  equal  to  the  sum  of  the  envelopes 
of  the  covariance  functions  for  each  mode.  The  maximum  of  the 
envelope  of  the  covariance  of  each  mode  occurs  at  a time  delay 
that  is  determined  by  the  source  direction,  detector  separation. 
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and  group  velocity  of  the  mode.  I suggest  that  It  say  be  nec- 
essary tc  average  the  square  of  starry  covariance  measurements 
to  Jeterclne  the  source  direction  with  time  delay  steering. 

Fven  so,  the  sun  of  the  different  sazlaa  of  the  envelopes  for 
the  different  directions  would  broaien  the  covariance  function. 
The  fluctuations  of  the  covariance  may  be  large  in  this  case 
and  correspondingly  the  apparent  noise  level  large. 

This  study  has  been  restricted  to  a noise  source  in 
a waveguide  with  infinite  layers.  The  layers  are  assumed  to 
have  constant  thickness  and  constant  properties.  If  the  layers 
are  not  uniform  in  the  area  of  the  array,  the  acoustical  pres- 
sures at  the  two  detectors  P1(t)  and  F2(t)  can  be  different. 
The  non- uniform  conditions  can  cause  the  covariance  function  to 
be  small  and  erratic.  Sder  these  conditions,  Measurement  of 
the  covariance  function  for  any  type  of  steering  would  indicate 
that  the  ~r -parent  noise  level  is  high. 

As  a final  consideration  let  us  assume  that  ve  have 
a simple  harmonic  source  in  the  waveguide.  The  frequency  of 
the  acoustical  pressure  is  the  same  for  all  detector  positions 
regardless  of  any  non-onifora  conditions  or  the  number  of  rodes 
propagating.  This  means  that  the  covariance  of  the  acoustical 
pressures  is  a function  of  the  phase  difference  of  tie  acousti- 
cal pressures  and  that  the  covariance  may  be  large,  whereas 
under  the  same  conditions  vith  a noise  scarce  having  small  band- 
width, the  covariance  would  tend  to  zero.  For  these  reasons  it 


is  difficult  to  predict  the  performance  of  a more  complex  array 
in  an  irregular  waveguide  on  the  basis  of  studies  made  with  a 
harmonic  source  and  a two-element  array. 
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